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SUMMARY 

The  theory  is  presented  of  a  technique  for  measuring  the  attitude  and 
position  of  cylindrical  objects  which  is  of  interest  in  studying  the 
behaviour  of  stores  ejected  from  aircraft  wing  pylons.  The  sight-line 
angles  of  the  tangents  to  the  top  and  bottom  edges  of  the  store  are  measured 
in  two  vertical  planes  using  linear  CCD  array  cameras,  and  equations  relat¬ 
ing  these  angles  to  the  required  attitude  and  position  are  derived  in  this 
Report.  The  measurement  of  bank  angle  from  the  observation  of  marks  painted 
on  the  side  of  the  store,  is  also  considered.  The  errors  in  position  and 
attitude  determination,  which  arise  from  uncertainties  in  the  angle  measure¬ 
ments  and  camera  geometry,  are  discussed  and  simple  first  order  approxima¬ 
tions  are  derived  for  these. 
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I  INTRODUCTION 

A  new  method  o£  measuring  the  position  and  attitude  of  cylindrical  stores  during 
their  release  or  ejection  from  aircraft  wing  pylons  is  currently  being  investigated  by 
Instrumentation  and  Trials  Department  RAE.  The  aim  is  to  make  available,  within  a  few 
minutes  of  store  release,  information  such  as  the  pitch  rate  and  separation  velocity 
time-histories,  thereby  eliminating  the  lengthy  time  scales  and  tedious  analysis 
associated  with  existing  photographic  techniques.  Initial  interest  is  focussed  on  the 
ground  testing  of  ejector  performance  but  the  wider  problems  of  airborne  carriage  and 
release  trials  are  being  kept  firmly  in  mind. 

The  technique  being  studied  employs  two  solid  state  cameras  incorporating  linear 
charge  coupled  imaging  devices  (CCDs)  to  measure  the  angles  of  the  sight-line  tangents  to 
the  top  and  bottom  edges  of  the  store,  in  two  separate  vertical  planes.  These  four 
measurements,  together  with  a  knowledge  of  the  store  radius  and  camera  positions,  enable 
the  position  of  the  store  in  two  dimensions  and  its  heading  and  elevation  attitude  to  be 
determined,  provided  that  the  cross-section  of  the  store  is  cylindrical  in  the  measurement 
zone.  It  follows  from  this  restriction  that  measurement  of  the  sight-line  tangent  angles 
is  incapable  of  determining  roll  about  the  longitudinal  axis  of  the  store,  or  motion 
along  this  direction.  In  principle  however,  bank  angle  may  be  determined  by  the  observa¬ 
tion  of  suitable  markings  painted  along  the  side  of  the  store  with  the  same  CCD  cameras 
as  used  for  measuring  heading  and  elevation  angles. 

The  aim  of  this  Report  is  to  present  the  geometrical  relationships  between  the 
store  position  and  attitude  and  the  measure  sight-line  tangent  angles  from  which  the 
former  may  be  calculated.  Although  not  required  in  the  present  application  the  computa¬ 
tion  of  bank  is  also  included  for  completeness.  Finally  the  influences  of  measurement 
and  alignment  errors  on  the  calculated  store  position  and  attitude  are  discussed.  The 
detailed  design  and  results  of  the  evaluation  of  an  experimental  system  will  be  reported 
elsewhere  in  due  course. 

An  advantage  of  the  technique  considered  here,  in  addition  to  its  fast  response 
time,  is  that  no  modification  of  the  store  is  necessary,  other  than  a  coat  of  paint  to 
provide  a  high  contrast  image  against  its  background.  A  major  disadvantage  however  is 
the  need  for  a  cylindrical  cross-section  at  least  in  the  regions  viewed  by  the  cameras. 

In  principle  this  could  be  overcome  by  the  use  of  additional  cameras  and  by  making  use  of 
the  known  profile  of  the  store,  however  the  computational  load  would  increase  considerably 
and  a  more  versatile  approach  in  this  instance  would  be  to  use  shape  correlation 
algorithms  on  a  complete  image  of  the  store.  This  would  probably  require  a  large  main¬ 
frame  computer  for  the  necessary  processing  speed  and  power.  In  contrast,  the  computation 
involved  in  the  technique  described  here  is  within  the  capabilities  of  a  hand-held 
programmable  calculator  and  since  at  least  two  of  the  stores  in  question  do  have 
cylindrical  cross-sections,  over  some  part  of  their  profile,  the  technique  clearly 
merits  further  consideration. 

The  term  attitude  is  used  in  a  restricted  sense  in  this  Report  to  include  only 
elevation  and  heading  angles.  The  third  parameter,  bank  angle,  normally  encompassed  in 
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the  definition  of  attitude  is  only  referred  to  in  the  specific  section  dealing  with  the 
derivation  of  the  bank  angle  equations  (section  2.5). 

2  THE  GEOMETRICAL  EQUATIONS 

2. 1  General 

The  basic  arrangement  of  two  linear  array  cameras  positioned  to  one  side  of  the 
store  whose  attitude  and  position  are  required,  is  illustrated  in  Fig  1.  Both  cameras 
are  accurately  aligned  so  that  each  CO)  array  has  a  fan  shaped  field  of  view  in  a 
vertical  plane  which  is  of  sufficient  extent  to  include  the  top  and  bottom  edges  of  the 
store  throughout  the  required  displacement  range.  If  the  store  is  cylindrical,  then  its 
cross-section  in  the  plane  of  either  camera  will  be  an  ellipse,  in  general,  whose  orienta¬ 
tion  and  eccentricity  are  functions  of  the  store's  attitude.  The  angles  of  the  sight¬ 
line  tangents  from  the  cameras  to  the  elliptical  sections  will  therefore  be  functions  of 
this  attitude  as  well  as  relative  position. 

The  two  cameras  are  assumed  to  lie  in  the  same  horizontal  plane  and  the  line 
joining  them  is  conveniently  taken  to  define  the  datum  heading  angle  relative  to  which 
store  heading  $  *  is  to  be  measured.  The  cameras  need  not,  of  course,  be  set  up  to 
view  exactly  at  right  angles  to  this  direction,  nor  need  their  viewing  planes  be  parallel 
to  each  other,  although  the  geometrical  relationships  will  obviously  simplify  in  these 
special  cases.  The  general  case  is  considered  here  however,  by  specifying  separate 
heading  angles  for  the  viewing  planes  of  the  two  cameras. 

In  the  analysis  which  follows,  the  equation  of  the  elliptical  section  of  the  store 
in  the  plane  of  one  of  the  cameras  is  derived  first  and  then  the  equations  of  the  sight¬ 
line  tangents  to  this  ellipse  from  the  camera  position  are  determined,  giving  the  sight¬ 
line  angles  as  functions  of  store  attitude  and  position.  Manipulation  of  the  sum  and 
differences  of  the  angles  to  the  top  and  bottom  edges  of  the  store  enables  a  set  of 
equations  to  be  derived  which  express  store  position  as  a  function  of  these  compound 
angles  and  attitude.  A  further  set  of  equations  is  derived  which  link  store  attitude 
to  store  position  and  the  known  geometry  of  the  cameras.  These  two  sets  of  equations 
must  then  be  solved  simultaneously  if  store  attitude  and  position  are  to  be  found 
independently.  Unfortunately  it  has  not  proved  possible  to  solve  these  equations 
analytically  but  an  interative  solution  is  presented  which  is  suitable  for  implementation 
on  a  confute  r  or  programmable  calculator  and  the  performance  of  this  technique  on  a  TI59 
calculator  is  discussed  in  section  3. 

The  equations  required  for  the  calculation  of  bank  angle  are  derived  in  section  2.5 
and  require  an  a  priori  knowledge  of  pitch,  yaw  and  position. 

2.2  Equation  of  the  store  section 

The  equation  of  the  store  section  in  the  camera  viewing  plane  can  be  found  from 
the  intersection  of  this  plane  with  the  three-dimensional  equation  of  the  cylindrical 
surface  of  the  store. 


*  For  a  complete  list  of  syrtols  see  page  34 
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If  the  store  has  a  radius  a  and  an  instantaneous  attitude  described  by  elevation 
angle  6  and  heading  ifi  ,  then,  for  the  spherical  polar  coordinate  system  shown  in 
Fig  2,  the  equation  of  its  three-dimensional  surface  is  shown  in  Appendix  A  to  be: 


-  jsin  a  sin  0  +  cos  a  cos  8  cos(0  - 


(1) 


where  R  -  radius  from  the  origin  of  a  point  P  on  the  store  surface 
a  -  elevation  angle  of  point  P  above  the  horizontal  plane 
8  -  bearing  of  point  P  in  the  horizontal  plane. 

Note  that  the  origin  of  the  coordinate  system  used  here  has  been  chosen  to  lie  on  the 
centre  line  of  the  store.  If  we  further  assume  that  the  viewing  plane  also  passes 
through  the  origin,  with  heading  angle  ,  then  the  store  section  in  this  plane  is 
given  simply  by  substituting  6  «  ^  in  equation  (1)  above.  We  may  also  rearrange  the 
trigonometric  functions  in  the  denominator  of  the  equation  to  give  finally  for  the 
equation  of  the  store  section: 

o  „2 

R2  -  - * -  (2) 

I  -  K  sin  (a  +  y) 


2  2 

where  K  *•  1  -  cos  0  sin  -  i(/)  ' 

y  -  tan  *|cot  0  cos  -  i|i)j  . 

Appendix  B  shows  that  this  represents  an  ellipse  with  minor  and  major  semi-axes 
of  a  and  <z/(l  -  K) ^  respectively,  rotated  through  an  angle  y  ,  where  y  is  the 
angle  between  the  minor  axis  of  the  ellipse  and  the  horizontal  plane. 

2.3  The  sight-line  tangents 

The  angles  of  the  sight-line  tangents  to  the  top  and  bottom  edges  of  the  store  may 
now  be  determined  by  finding  the  equations  of  the  tangents  to  the  elliptical  store 
section  described  by  equation  (2)  above.  This  operation  is  facilitated  by  changing  the 
origin  of  the  coordinate  system  employed  from  the  centre  of  the  store  to  the  camera 
position,  whereupon  the  tangents  become  simple  radial  vectors  passing  through  the  new 
origin. 

If  the  position  of  the  store  centre  is  specified  by  range  Rq  and  depression 
angle  aQ  in  the  new  coordinate  system  as  shown  in  Fig  3  then  the  new  coordinates  of  a 
point  P  on  the  surface  of  the  store  R’ ,  a'  may  be  related  to  the  old  coordinates 
R,  a  by  the  equations: 


R  cos  a  *  R'  cos  a'  - 


cos  a„ 


R  sin  a 


Rq  sin  <Jq  -  R'  sin  a' 


(3) 


The  angles  a'  and  <Sq  have  been  defined  as  depression  angles  here  in  anticipation  of 
the  likely  experimental  configuration. 


m 
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Transforming  the  store  section  into  the  new  coordinate  system  using  equation  (3) 
gives,  after  some  manipulation: 

R,2{l  ~  K  sin2(a'  -  y)^  -  2R'R0|cos(a0  -  a')  -  K  sin(aQ  -  y)  sin(a'  -  y)| 

+  R2{t  -  K  sin2(a0  -  y)|  -  a2  -  0  .  (4) 

This  is  a  quadratic  in  R'  indicating  that  for  a  given  radius  vector  at  depression 
angle  o'  there  are  in  general  two  points  of  intersection  with  the  elliptical  store 
section.  As  the  vector  approaches  a  tangent  position  however  the  two  roots  of  the 
equation  must  converge  to  a  single  value  and  this  criterion  of  equal  roots  can  be  used 
to  establish  the  angles  of  the  tangents. 

2 

The  roots  of  a  quadratic  of  the  general  form  Ax  +  Bx  +  C  ■  0  are  equal  when 

2 

B  »  4AC  ,  hence  from  equation  (4)  above  the  tangent  angles  must  satisfy: 

Ro{COS(“o  ~  “')  “  K  ain(aQ  ~  y)  sin(a'  -  y)^ 

-  jl  ~  K  sin2(a'  -  Y)}{&o(l  -  K  sin2(aQ  -  y))  -  a2}  .  (5) 

In  order  to  solve  this  for  the  angle  a'  we  can  multiply  out  and  make  the  substitutions: 
cos -  o')  -  cos j(aQ  -  y)  -  (o'  -  y)^ 

-  cos(aQ  -  y)  cos (a'  -  y)  +  sin(aQ  -  y)  sin(a'  -  y) 

cos2 (a'  -  y)  «  j(l  +  cos  2(u'  -  y)J 
sin2(a'  -  y)  ■  J  (l  -  cos  2(a*  -  y)) 

and 

cos (o'  -  y)  sin(a'  -  y)  ■  J  sin  2(a'  -  y) 
which  yields  the  equation: 

P  cos  2 (o'  -  y)  +  Q  sin  2(a'  -  y)  -  S  (6) 

2 

whore  P  -  (I  -  K)  cos  2(aQ  -  y)  + 

*6 

Q  -  (1  -  K)  sin  2(aQ  -  y) 

_2 

S  -  1  -  K  -  (2  -  K)  . 

*0 
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However,  this  is  recognised  to  be  of  a  standard  trigonometric  form  enabling  us  to  write: 


(p2  +  Q2)  sin|2(a'  -  y)  +  5}  -  S 


(7) 


where  £ 


Rearranging,  we  find  that  the  angle  of  depression  of  the  tangent  is  given  by: 


a' 


(8) 


The  two  tangent  angles  required  are  given  by  the  two  values  of  the  sin  *  term 
between  0  and  2ir  which  are,  of  course,  simply  related  since  their  sum  is  equal  to  ir 
It  follows  therefore,  that  if  the  tangent  angles  to  the  top  and  bottom  of  the  store  are 
a j  and  a j  respectively,  then: 


and 


(9) 


(10) 


These  sum  and  difference  equations  may  be  further  developed,  as  shown  in  Appendix  C 
to  give  the  following  equations  for  the  range  and  depression  angle  of  the  store  centre: 


(al  + 


o2)  -  sin-'jy^^  sin{2y  -  (a,  ♦  a2>}J 


01) 


2  ♦  k|cos2(o2  -  otj)  cos  2(<*q  -  y)  -  -  2  cos(a2  -  Oj)  x 

*  | I  -  K  sin2(Og  -  y)^2  -  K  +  R  cos2(oQ  -  y)  cos2(a2  -  »,)}| 
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The  ranges  and  depression  angles  of  Che  store  centres  in  the  two  camera  planes  will 
of  course  be  different  in  general  but  can  be  related  for  a  given  store  attitude  by 
considering  the  known  camera  geometry. 

2.4  Camera  geometry 

The  two  sets  of  parameters  which  apply  to  the  two  different  camera  planes  may  be 
conveniently  distinguished  by  the  use  of  primed  and  unprimed  symbols.  Thus  the  set  of 
unprimed  synbols  used  up  to  now  {rq,Oq,Oj  »a2»^o  ^  be  taken  as  applying  to  the  forward 
camera,  while  the  corresponding  primed  symbols  ^o,aO*a'l  ,a2*^o}  re^er  to  tlie  camera. 

Appendix  0  shows  how  the  range  and  depression  angles  in  any  two  camera  planes  are 
related  for  a  given  store  attitude  (6,t)0  and  camera  separation  D  . 


Thus: 


tan  (*q 


lain  6  cos  aQ  sin(4i^ 


,)  +  cos  0  sin  a„  sin(i)j’  -  <10  >  +  D  sin  0  sin  <|)i 


cos  ©{Rq  cos  aQ  sin(ij*0  -  ♦)  *■  D  sin  <|i| 


R_  cos  a_  sin(<Ji_  -  <|0  -  D  sin 

"6  -  - cos  a^sfn??»— 3 -  *  (,4> 

The  term  in  may  be  eliminated  from  equation  (14)  by  substitution  from  equation  (13) 

but  this  is  tedious  and  gives  a  very  cumbersome  result,  whereas  the  above  form  is 
perfectly  adequate  for  numerical  evaluation,  provided  that  equation  (13)  is  evaluated 
first. 

Appendix  D  also  shows  that  attitude  may  be  determined  from  a  knowledge  of  the  range 
and  depression  angle  of  the  store  centre  in  the  two  camera  planes: 


*0  cos  o0  sin  4>0  ~  Rq 
D  +  Rq  cos  aQ  cos  d>Q  “  i 


-  %  g°8  ao  8llt  *0 

*0  ’  COS  Og  cos 


tan  0  *  sin 


%  8ip  ao  -  *Q  8in  a0 

«0  cos  a0  8in  *0  "  *b  cos  ao  sin  *0 


The  term  in  ^  in  equation  (16)  may  be  eliminated  by  substitution  from  equation  (15)  but 
again  the  above  simpler  form  is  adequate  for  numerical  evaluation,  provided  that  41  is 
determined  first.  The  second  pair  of  equations  is  of  course  equivalent  to  the  first 
(equations  (13)  and  (14)),  having  been  derived  from  them  merely  by  a  rearrangement  of 
terms,  and  collates  the  set  of  equations  necessary  for  the  determination  of  store 
position  and  attitude  in  elevation  and  heading.  The  solution  of  the  equations  is 
discussed  in  section  3. 

It  should  be  noted  that  equations  (13)  to  (16)  inclusive  can  be  divided  top  and 
bottom  by  store  radius  leaving  only  normalized  dimensions  Rq/<z  and  D/a  and  showing 
that  behaviour  can  be  simply  scaled  according  to  store  radius. 
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2.5  Bank  angle 

Although  not  required  for  the  present  application  the  measurement  of  bank  angle  is 
presented  here  for  completeness.  Its  determination  requires  an  a  priori  knowledge  of 
the  store  position  and  attitude  in  heading  and  elevation,  whereupon  it  may  be  directly 
calculated  from  the  observed  sight-line  angles  to  marks  painted  on  the  side  of  the  store. 

It  is  shown  in  Appendix  E  that  the  sight-line  depression  angle  ,  in  the  plane  of 
the  forward  camera,  of  longitudinal  marks  on  the  side  of  the  store  at  a  bank  angle  $  is 
given  by: 


*0  . 


tan  a„ 


—  sin  <*Q  cos  0  sin(ipQ  -  ip)  -  sin  $  sia(ipQ  -  <p)  +  sin  6  cos  $  cos(<»q  -  <») 


cos  0 


—  cos  ctg  sin(fyj  -  i/i)  -  cos 


(17) 


This  may  be  rearranged  to  give  bank  angle  <p  as  a  function  of  the  measured  angle  : 


*  ” 


•  -1, 

n  -  sin  < 


*0  /  1 
—  cos  0  sin(ij)Q  -  tfi)<cos  tan  -  sin  oQv 


+  cos20|tan2a3  +  2  tan  tan  0  cos(s pQ  -  ip)  -  cos2(ipQ  -  V>)|J 


OB) 


j(cos  0  tan  +  sin  0  cos(i|/q 
where  r,  -  tan"  \ - kln<^'-  - 


Similar  expressions  apply  to  the  aft  camera  of  course  but  with  Rq,  aQ  and  ^ 
replaced  by  R^,  o^  and  ip^  respectively.  The  direct  calculation  of  a  mean  roll  angle 
in  this  way,  is  quire  straightforward  once  store  position  and  attitude  in  elevation  and 
heading  are  known. 

3  SOLVING  FOR  POSITION  AND  ATTITUDE 

Evaluating  equations  (11)  and  (12)  in  the  two  camera  planes  (using  primed  and 
unprimed  symbols  as  appropriate),  together  with  a  pair  of  the  camera  geometry  equations 
from  section  2.4,  provides  a  set  of  six  equations  altogether,  giving  the  parameters  Rq, 
V  “o*  °o»  9  and  ♦  in  terms  of  each  other  and  known  or  measured  quantities,  and  which 
in  principle  can  be  solved  to  find  these  parameters.  An  analytical  solution  of  these 
equations  by  the  conventional  technique  of  substitution  and  elimination  would,  however, 
be  extremely  tedious  even  if  it  were  possible  and  this  route  has  not  been  seriously 
considered.  Instead  a  relatively  si^>le  iterative  technique  has  been  employed  which  is 
easily  implemented  on  a  programmable  calculator  or  microcomputer  and  since  the 
computational  load  is  sufficient  to  rec*  re  the  •"  of  such  a  machine  in  any  case,  this 
does  not  present  any  additional  problems.  ’  ju  *  erative  procedure  adopted  is  described 
below. 
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The  £irst  step  in  the  procedure  is  to  make  an  initial  estimate  of  the  ranges  and 
depression  angles  of  the  store  centre  in  the  two  camera  planes.  This  is  conveniently 
achieved  by  assuming  that  the  store  section  in  these  planes  is  circular.  that  there 
is  no  pitch  (8  -  0)  and  the  camera  planes  are  orthogonal  to  the  store  axis 
((^q  -  4>)  -  (i|£  -  iji)  «  ir/2).  In  this  special  case  equations  (11)  and  (12)  simplify 
considerably  giving,  for  the  forward  camera: 


ao 

i(«,  +  o‘2) 

(19) 

h  . 

(**2  ~  a\\ 
cosecl - = - 1 

(20) 

a 

\  2  / 

and  sioiilarly  for  the  aft  camera: 

* 

(21) 

K  , 

T2  ‘  “!\ 

cosecl - s - 1 

(22) 

a 

\  2  / 

Using  these  initial  estimates  for  the  positions  of  the  store  centres  in  the  two 
camera  planes  it  is  now  possible  to  obtain  a  more  realistic  estimate  of  the  store 
attitude  knowing  the  camera  geometry  and  using  equations  (15)  and  (16)  which  are 
repeated  here  in  their  non-dimensional  form. 


R, 


0 

T  C08  ao  8itl  *0  "  T  cos  ao  sin  *0 

tan  i|«  -  - g- - gi - 

£  +  T  cos  aQ  cos  *0  -  cos  cos 


(23) 


tan  8  ■  sin 


R0  .  ,  ^  . 

—  8in  “o  '  T  sin  ao 

rT 


cos  aQ  sin  <|»0  -  —  cos  sin 


(24) 


This  new  estimate  of  attitude  may  in  turn  be  used  to  obtain  a  better  estimate  of 
the  true  position  of  the  store  centres  in  each  of  the  camera  planes.  Thus  from 
equations  (11)  and  (12)  for  the  forward  camera: 


•  i[(al  +  a2>  “  *in~‘  sin{2r  -  (01  +  a2)}[] 


(25) 


£ 

4 


2  +  k|cos2 * 4(o2“  Oj)  cos  2(«q  -  y)  -  l|  -  2  cos(a2  -  Oj)  x 

x  |l  -  K  sin2 (0^  -  y)^2  -  K  +  K  cos2(n0  -  y)  cos2(a2  -  aj)|j 


4  *  j  R -R  sin2(o2  -  Oj) 


..  (26) 


where  K 


1  i 


2  2 

1  -  cos  8  sin  (ipQ  -  iji) 


Y  ■  tan  '  |cot  8  cos(ipQ  -  iji)^ 
and  similarly  for  the  aft  camera 


[<*i  +  a2)  "  sin-'j-p^.^  sin{2Y'  -  (a;  ♦ 


(27) 


2  +  K'|cos2(a2  “  a'|)  cos  2(a^  -  y')  -  1 1  “  2  cos(a2  ~  aP  X 

*  jl  -  K  sin2(a^  -  Y')|2  -  K'  +  K'  cos2(a^  -  y')  cos2(a^  -  a j) J 


2 

4  +  j~  Ki-  sin2  (a^  -  aj) 


(28) 


2  2 

where  K'  -  1  -  cos  8  sin  (iji^  -  ip) 

Y*  ■  tan  *|cot  8  cos(ijjg  -  i|()^ 

This  procedure  may  be  repeated  by  the  cyclic  evaluation  of  equations  (23)  to  (28) 
until  the  desired  accuracy  is  achieved. 

A  programme  to  implement  the  procedure  described  above  has  been  written  for  a  TI59 
programmable  calculator  and  has  an  execution  time  of  approximately  1  min/iteration.  It 
is  estimated  that  this  same  programme  would  run  some  7  times  faster  on  a  small  desk  top 
machine  such  as  the  Hewlett  Packard  KP85. 

The  number  of  iterations  required  before  the  answer  is  obtained  to  a  given  degree 
of  accuracy  depends  on  the  store  attitude  and  position  and  generally  increases  markedly 
with  the  eccentricity  of  the  store  section  in  the  camera  planes.  It  also  depends  on  the 
orientation  of  this  elliptical  section  and  can  require  more  than  twice  the  number  of 
iterations  when  the  major  axis  of  the  ellipse  is  nominally  at  right  angles  to  the  mean 
camera  sight-line  compared  to  when  the  major  axis  is  aligned  with  this  direction. 

The  convergence  of  the  iterative  solution  is  illustrated  in  Tables  1  and  2  for  a 
particular  case  being  considered  where  the  cameras  are  separated  by  only  1.25  store  radii 
(D/a  ~  1.25)  and  view  at  90°  to  the  datum  heading  direction.  The  use  of  smaller  values 
of  D/a  than  this  is  unlikely  in  practice  and  to  this  extent  this  example  represents 
a  worst  case.  The  store  is  assumed  to  be  nominally  at  a  range  of  five  store  radii  and  at 
a  depression  angle  of  30°  but  with  elevation  and  heading  angles  of  10°  in  the  first  case 
and  30°  in  the  second.  It  is  seen  that  the  initial  estimates  of  attitude  (iteration 
>  number  0  in  the  tables)  are  within  0.5°  when  the  elevation  and  heading  angles  are  10°  but 

_  J 
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are  in  error  by  up  to  7°  for  the  30°  attitude.  Convergence  is  also  more  rapid  in  the  10° 
case  with  the  answer  being  correct  to  0.001°  in  only  two  iterations  while  ten  iterations 
are  required  in  the  30°  case.  The  data  given  in  the  tables  for  the  infinite  iteration 
case,  are  in  fact  the  initial  data  used  to  calculate  the  sight-line  tangent  angles  fed 
into  the  iterative  prograame  and  which  were  computed  using  equations  presented  in  this 
Report. 

The  convergence  behaviour  for  elevation  and  heading  angles  both  of  opposite  sign  to 
those  just  considered  would  be  identical  to  that  shown  above,  since  the  eccentricity  and 
orientation  of  the  store  sections  would  be  the  same,  but  convergence  would  be  faster  if 
only  one  were  of  different  sign,  since  the  ellipse  would  now  be  oriented  more  towards  the 
cameras  as  discussed  earlier.  Thus  for  elevation  8  =*  -30°  and  heading  i|»  =  30°  the 
initial  estimates  of  attitude  are  within  1.7°  and  the  solution  is  correct  to  0.001°  in 
only  four  iterations. 

4  ERROR  ANALYSIS 

4. 1  General 

The  foregoing  equations  may  be  used  to  compute  the  detailed  effects  of  measurement 
errors  occurring  in  the  sight-line  tangent  angles  (a ^  ,°2,aj  ,ap  or  errors  in  setting 
up  the  camera  geometry  (*j'g,i|'g,D)  for  any  particular  store  position  and  attitude.  However 
it  is  more  useful  to  the  system  designer  to  have  simple  analytical  expressions  which 
relate  the  resultant  errors  in  position  and  attitude  to  the  measurement  or  setting  up 
errors  where  these  can  be  derived. 

An  estimate  of  the  effect  of  finite  or  limited  camera  resolution  is  of  particular 
interest  here.  In  practice  resolution  may  be  determined  by  the  size  of  the  individual 
CCD  photosites  and  the  focal  length  of  the  lens  employed,  giving  rise  to  a  quantization 
in  the  angular  measurements  and  it  is  important  to  know  how  the  errors  associated  with 
this  translate  into  uncertainties  in  store  attitude  and  position.  The  effect  of 
measurement  errors  has  therefore  been  considered  in  some  detail.  Errors  in  the  determina¬ 
tion  of  bank  angle  have  not  been  considered. 

4.2  Measurement  errors 

Simple  analytical  expressions  for  the  effect  of  measurement  errors  may  be  found  by 
restricting  attention  to  the  moderate  pitch  and  yaw  situation,  eg  such  that 
sin  0  *  tan  0  *  0  ,  to  when  the  store  is  some  distance  away  from  the  cameras 
((Rg/a)2  ►  1 )  ,  and  to  when  the  cameras  are  orthogonal  to  the  reference  heading  direction 
(*0  “  ♦q  ■  */2)  .  To  a  first  approximation  the  effects  of  measurement  errors  can  then  be 
assumed  to  be  independent  of  actual  store  attitude  and,  more  specifically,  these  effects 
can  be  quantified  by  considering  the  zero  pitch  and  yaw  case,  which  greatly  simplifies 
the  mathematics  involved. 

The  uncertainties  in  range  6Rg  and  depression  angle  Sag  can  be  related  to  the 
errors  5a j  and  602  in  the  measurement  of  the  sight-line  tangent  angles  a j  and  02 
respectively,  for  the  forward  camera,  by  the  partial  differential  equations  for  separable 
variables: 


13 


6R„ 


5a„  - 


3Ro  *  *V 

1^6al  3a^  6  Z 


3a, 


3a, 


'0 


3a,  1  3a2 


(29) 


(30) 


Similar  expressions  apply  to  the  aft  camera  using  the  appropriate  primed  variables. 

The  partial  derivatives  in  the  above  equations  may  be  found  in  the  zero  pitch  and 
yaw  case  (0  *  0,  ip  -  0)  by  differentiation  of  equations  (19)  and  (20),  substitution  of 
these  yields: 


and 


5  a, 


R0 

2§(5a2-6a.) 


J(6a,  +  Sa^) 


(31) 


(32) 


These  equations  together  with  their  counterparts  for  the  aft  camera  enable  the  mean 
positional  errors  in  range  and  depression  angle,  ARq  and  Aa^  respectively,  to  be 
calculated,  using: 

ARq  -  i(6R0  +  SRq)  (33) 


Aa, 


}(<$a0  +  S(Xq)  . 


(34) 


The  errors  6i|<  and  60  in  the  heading  and  elevation  angles  respectively,  may  be 
determined  similarly,  from  their  partial  differential  equations: 


and 


6*  -  6R  +  6Rl  +  3^-  6a  +  6a' 

3Rq  0  3Rg  V  3aQ  0  3aQ  0 


66  ’  3^  6R0  +  4  6R°  +  ^0  &a°  +  4  S  ‘ 


(35) 


(36) 


The  partial  derivatives  required  here  may  be  obtained  from  the  attitude  equations  (15) 
and  (16),  which  simplify  in  the  case  being  considered,  to: 


*  tan  i)<  * 


R0 


cos  a. 


% 


cos  a; 


(37) 


0  *  tan  0  * 


Rq  8in  a0  "  *0  8in  ac 


(38) 


Thus  for  Oq  **  Oq  and  R^ 


we  have: 


cos  a. 


6i(i  “ 


-  6Rq)  +  -ysin  (5°'  -  5<*q) 


(39) 
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se 


sin  a 

— r9-  (6Ro  “  6Ro> 


*0 


D  cos  oQ  (5«q  -  «<*0) 


(40) 


The  positional  errors  given  by  equations  (31)  and  (32)  and  their  counterparts  for  the  aft 
camera  may  be  substituted  in  the  above  to  give  the  attitude  errors  directly  in  terms  of 
the  measurement  errors  dOj,  6a'j  and  60^  • 

Of  particular  interest  to  the  modelling  of  the  effects  of  finite  camera  resolution, 
is  the  case  where  all  four  measurements  can  be  in  error  by  the  same  amount,  of  magnitude 
Aa  say,  but  where  the  sign  of  the  error  can  be  positive  or  negative  in  an  arbitrary 
manner.  Since  there  are  four  measurements  there  will  be  16  different  possible  error 
combinations  altogether,  but  only  eight  of  these  will  lead  to  different  magnitudes  of 
error  in  store  position  and  attitude,  the  other  eight  giving  similar  magnitudes  but 
opposite  signs.  It  is  a  relatively  simple,  if  somewhat  tedious  matter  to  evaluate  these 
eight  different  error  conbinations  but  only  those  giving  rise  to  the  worst  position  and 
attitude  errors  will  be  discussed  further  here. 


The  largest  range  error  occurs  when  both  measurements  to  the  top  of  the  store 
(<Xj  and  a|)  are  too  large  and  those  to  the  bottom  (a^  and  ap  are  too  small,  or 
vice  versa,  giving; 


A80 


(41) 


The  worst  case  depression  angle  error  obviously  arises  when  all  the  measurement  errors 
have  the  same  sign  and: 

AOq  »  ±  Aa  .  (42) 


The  largest  attitude  errors  generally  occur  when  the  measurement  errors  combine  to 
make  the  nose  of  the  store  appear  closer  to  the  camera  than  it  really  is  and  the  tail 
appear  further  away,  or  vice  versa.  Thus,  for  the  forward  camera,  the  top  sight-line 
angle  might  be  too  small  and  that  to  the  bottom  too  large,  while  for  the  aft  camera  the 
top  angle  would  be  too  large  and  the  bottom  one  too  small.  The  errors  in  this  situation 
are  given  by: 


6ip 


cos  Aa 
Da 


(43) 


and 


66 


2Rq  sin  o0 
Da 


Aa 


(44) 


For  purely  horizontal  viewing  (a^  -  0)  the  error  in  yaw  reaches  a  maximum  but  the 
pitch  error  66  is  then  zero  for  this  particular  error  format.  The  worst  case  pitch 
error  in  this  situation  occurs  when  both  measurement  errors  for  the  forward  camera  are  of 
the  seamt  polarity,  eg  making  the  nose  appear  high,  while  those  at  the  tail  are  both 
opposite  to  this,  making  the  tail  appear  low.  The  pitch  error  is  then  given  by: 

21L.  cos  an  Aa 

68  -  — S - fl—E -  .  (45) 
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It  is  easily  seen  that  this  exceeds  the  error  given  by  equation  (44)  above  for  small 
angles,  such  that  Oq  <  tan  '(a/ Rg)  ,  and  this  very  roughly  equates  to  the  store  centre 
being  less  than  one  store  radius  below  the  horizontal  plane  containing  the  camera,  ie  to 
small  drop  heights. 

Consideration  of  all  16  possible  measurement  error  combinations  gives  average 
position  and  attitude  errors  which  are  less  than  half  of  the  worst  case  values  presented 
above.  In  practice,  with  a  random  error  distribution,  there  will  be  occasions  when  some 
or  all  of  the  measurement  errors  are  zero  and  the  resultant  attitude  and  position  errors 
will  therefore  be  even  smaller.  It  is  common  practice  in  such  a  situation  to  equate  the 
worst  case  errors  with  three  times  the  standard  deviation  of  the  random  errors. 

The  worst  case  attitude  errors  can  be  many  times  greater  than  the  measurement 
errors  in  the  sight-line  tangent  angles  as  numerical  evaluation  of  equations  (43)  and 
(44)  quickly  shows.  Thus  in  the  specific  case  considered  in  section  3  (Rg/a  «  5, 

D /a  -  1.25,  ag  *  30°)  we  have : 

6 ip  "  34.6  4a 

(46) 

6e  -  20.0  Aa  . 

Hence  to  measure  attitude  to  ±1°  in  this  case  it  would  be  necessary  to  measure  the  sight¬ 
line  tangent  angles  with  an  error  of  less  than  ±0.03°  and  it  is  clear  that  very  stringent 
demands  can  be  placed  on  the  system  designer  in  this  way. 

The  major  questions  posed  by  the  system  designer  are  likely  to  involve  the 
optimisation  of  the  viewing  geometry,  eg  camera  spacing  and  range,  in  order  to  minimise 
the  influence  of  the  measurement  errors  discussed  above.  The  dependencies  contained  in 
the  error  equations  indicate  that  the  technique  is  best  suited  to  long  stores,  enabling 
the  cameras  to  be  placed  a  large  distance  apart  (D  large)  and,  at  first  sight,  suggest 
greatly  improved  performance  at  relatively  short  ranges  (Rg  small). 

Clearly  camera  separation  must  be  made  as  large  as  practically  feasible  although 
it  is  really  the  separation  of  the  viewing  planes  at  the  store  position  which  is  of 
inq>ortance  here  and  so  toe-in  or  out  of  the  cameras  can  be  used  as  appropriate,  with 
benefit,  when  the  choice  of  camera  position  is  otherwise  constrained. 

The  apparent  dependence  of  attitude  errors  on  the  square  of  range  in  equations  (43) 
and  (44)  can  be  substantially  modified  by  the  interdependencies  of  depression  angle  Og  , 
angular  resolution  Aa  and  range  Rg  which  arise  in  practice.  The  designer  ie  usually 
interested  in  observing  the  motion  of  the  store  as  it  falls  through  a  predetermined 
vertical  drop  height  H  and  obviously  as  this  is  viewed  from  greater  and  greater  range 
the  sight-line  depression  angles  involved  become  smaller  and  hence  ag  which  enters 
into  the  error  equations  is  also  smaller.  Furthermore  the  focal  length  of  the  camera 
lens  can  be  increased  as  it  moves  away  to  maintain  the  field  of  view  in  the  plane  of  the 
store  equal  to  the  required  drop  height.  This  also  maintains  a  constant  spatial 
resolution  in  the  store  plane,  at  least  to  a  first  approximation,  if  it  is  assumed  that 


resolution  is  limited  by  the  CCD  element  size,  but  angular  resolution,  Aa  ,  till  then 
actually  i strove  in  inverse  proportion  to  the  increase  in  range  Rq  . 

The  result  of  such  interdependencies  is  that  the  error  in  heading  angle,  6i|i  ,  now 
varies  in  direct  proportion  to  the  horizontal  separation  of  the  camera  and  store 
(Rq  cos  Oq)  rather  than  with  Rq  ,  while  the  errors  in  elevation  angle  66  ,  given  by 
equation  (44)  for  large  (Xq  ,  become  directly  proportional  to  the  vertical  separation  of 
the  camera  and  store  (Rq  sin  c»q) .  The  latter  can  be  minimised  by  arranging  for  near 
horizontal  viewing  whenever  possible  and  in  this  event  the  elevation  errors  depend  only 
on  the  vertical  drop  height  H  which  must  be  covered.  As  indicated  earlier  this  will 
usually  be  a  fixed  parameter  in  the  experiments  and  the  worst  case  elevation  angle  error 
then  becomes  essentially  independent  of  range. 

In  the  small  angle  case  when  equation  (43)  applies  (H  <  a)  elevation  errors  will 
increase  with  range  due  to  the  associated  decrease  in  <Xq  .  This  variation  will 
typically  be  small  since  c»q  is  small  anyway. 

Different  dependencies  in  range  will  of  course  be  obtained  if  the  designer  invokes 
different  practical  links  between  resolution,  range  and  depression  angle  to  those  assumed 
above.  It  is  left  to  the  reader  to  evaluate  these  as  necessary.  The  general  conclusions 
to  be  reached  here  however,  are  that  the  influence  of  measurement  errors  on  attitude 
can  be  minimised  by  wide  camera  separation,  near  horizontal  viewing  and,  in  the  typical 
case  considered  (Aa  «  1/Rq),  minimising  viewing  range,  although  the  benefits 
attributable  to  the  latter  constraint  are  less  than  would  appear  at  first  sight  to  be  the 
case. 

Changing  the  camera  separation  D  or  nominal  viewing  angle  Oq  is  seen  to  have 
no  effect  on  the  magnitude  of  the  worst  cases  position  errors  (ARq  and  AOq)  arising 
from  measurement  errors  in  the  sight-line  tangent  angles,  as  might  be  expected.  In 
discussing  the  influence  of  viewing  range  on  these  position  errors  however,  we  may  again 
invoke  practical  relationships  between  range  depression  angle  and  angular  resolution, 
thus  modifying  the  apparent  dependencies  given  in  equations  (41)  and  (42).  Thus  if  we 
again  vary  lens  focal  length  to  give  an  angular  resolution  An  inversely  proportional  to 
range,  we  find  that  the  error  in  range  ARq  is  now  directly  proportional  to  the  range 
Rq  ,  while  the  angular  error  AoQ  varies  inversely  with  this.  The  net  effect  is  that 
positional  uncertainty  along  the  line  of  sight  is  proportional  to  range  while 
uncertainty  at  right  angles  to  the  sight-line,  given  by  the  product  RqAOq  ,  is 
essentially  constant.  For  near  horizontal  viewing  conditions  these  directions  may  be 
roughly  equated  with  the  local  horizontal  and  vertical  respectively  so  that  in  this  case 
horizontal  position  error  will  increase  in  direct  proportion  to  viewing  range,  while  the 
uncertainty  in  height  remains  constant.  Once  again  therefore  the  overall  balance  is 
probably  in  favour  of  shorter  viewing  ranges. 

The  influence  of  store  radius  on  the  effects  of  measurement  errors  is  seen  in 
equations  (41),  (43)  and  (44)  which  indicate  a  reduction  in  range  and  attitude  errors 
for  large  radius  stores,  all  other  things  being  equal.  Physically  this  goes  together 
with  the  requirement  for  shorter  viewing  ranges  noted  above  and  the  ability,  observed 
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in  section  2,  to  normalise  all  the  dimensions  involved  in  the  geometrical  equations  by 
store  radius,  which  in  themselves  would  suggest  that  normalised  range  R^/a  is  the 
important  parameter  here  and  the  one  to  be  minimised. 

4.3  Setting  up  errors 

The  effects  of  errors  in  setting  up  the  camera  geometry,  both  in  attitude  and 
separation,  may  be  discussed  qualitatively  in  general,  or  confuted  accurately  for 
specific  cases  using  the  iterative  technique  described  earlier. 

Errors  in  the  camera  separation  D  would  not  be  expected  to  lead  to  significant 
errors  in  store  position,  at  least  to  a  first  order  approximation  when  near  orthogonal 
viewing  conditions  are  used.  However  errors  in  store  attitude  can  be  generated  and  these 
would  be  expected  to  vary  in  direct  proportion  to  the  errors  in  D  .  Since  the  latter 
are  easily  limited  in  typical  practical  situations  to  much  less  than  1Z,  this  error 
source  can  be  effectively  eliminated. 

The  effects  of  errors  in  the  camera  heading  angles  (<Pq  and  i^)  can,  however,  be 
much  more  serious,  especially  when  the  camera  separation  is  small  and  the  viewing  range 
large,  since,  as  indicated  earlier,  it  is  really  the  separation  between  the  camera 
planes  at  the  store  which  is  important  in  determining  attitude.  It  is  easy  to  deduce  an 
expression  for  the  fractional  change  in  store  attitude  resulting  from  a  change  in 

both  camera  heading  angles,  the  worst  case  being  when  both  cameras  angle  in  towards 
each  other  or  both  toe  out.  To  a  first  approximation  the  attitude  errors  would  then  be 
expected  to  be  given  by: 

<56  <54*  2R0CO8°0,,  ,,,, 

T  "  T  - 5 - 6,"o  *  (47) 

These  errors  can  obviously  be  minimised  by  making  the  camera  separation  D  as 
large  as  possible  and  selecting  a  small  horizontal  working  distance  (Rq  cos  .  In  the 
particular  case  considered  earlier,  where  Rg/a  ■  5,  D/a  »  1.2S  and  -  30°  , 
equation  (47)  shows  that  to  measure  attitude  to  ±1°  over  a  range  of  ±20°  in  store 
elevation  and  heading,  camera  heading  angle  must  be  known  to  better  than  t0.4°. 

Camera  heading  angle  errors  would  be  expected  to  have  a  negligible  effect  on  store 
position,  to  a  first  approximation,  for  near  orthogonal  viewing  conditions 
((<l»0  -  *)  *  (1>q  -  ♦)  *  w/2)  . 

Errors  in  the  camera  elevation  angles  will  lead  to  errors  in  the  determination  of 
the  store  depression  angle  Oq  if  these  are  in  the  same  direction  for  both  cameras  and 
to  errors  in  store  attitude  if  of  opposite  polarity.  The  principal  attitude  error  would 
be  in  elevation  for  near  horizontal  viewing.  There  is  of  course  a  close  parallel  here 
with  some  of  the  measurement  error  combinations  discussed  in  section  4.2  and  similar 
equations  are  obtained  for  modelling  purposes.  Thus,  by  analogy  with  equations  (42) 
and  (45),  the  worst  case  depression  angle  error  %  resulting  from  a  camera  elevation 
angle  error  £9  is: 

Aan  -  ±  £0 
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which  for  near  horizontal  viewing  trill  mainly  result  in  an  error  in  store  height 
measurement,  and  the  worst  case  store  elevation  angle  error  £8  is  given  by: 


£6 


28°  “•  *2  se 


The  associated  store  heading  angle  error  is: 


Sip 


2R0  8in  °0 
D 


£0 


which  tends  to  zero  for  horizontal  viewing  (a^  23  0) . 

The  resultant  errors  in  attitude  and  height  measurement  will  all  increase  with 
greater  viewing  range  and  the  attitude  errors  will  also  become  worse  for  smaller  camera 
separations.  For  R^/a  -  5,  D/a  »  1.25  and  <*q  -  30°  the  elevation  and  heading  errors 
become  respectively  seven  and  four  times  the  error  in  camera  elevation  angle  and  hence 
for  attitude  to  be  determined  to  ±1°  camera  elevation  angle  must  be  accurate  to  ±0.14°. 

Finally,  uncertainties  in  the  radius  of  the  store  would  be  expected  to  result 
mainly  in  proportional  errors  in  the  store  range  with  a  minimal  effect  on  attitude,  any 
taper  in  the  store  profile  between  the  two  camera  planes  would,  however,  be  much  more 
serious  in  this  respect. 

4.4  Discussion 

It  is  clear  from  the  error  equations  presented  in  sections  4.2  and  4.3  that  the 
effects  of  both  measurement  errors  and  setting  up  errors  can  be  minimised  in  general  by 
working  at  relatively  close  range  and  with  as  large  a  separation  of  the  camera  planes  as 
possible.  While  the  effects  of  measurement  errors  will  usually  be  random  in  nature, 
setting  up  errors  will  be  systematic  and  are  therefore  less  of  a  problem  in  practice 
since  they  can  be  largely  eliminated  by  calibration  and  they  do  not  in  any  case 
significantly  affect  the  measurement  of  store  velocities  or  angular  rates  which  are 
frequently  the  variables  of  prime  interest. 

Calculations  of  the  errors  in  specific  cases,  using  the  iterative  techniques 
described  in  section  3,  has  confirmed  the  general  validity  of  the  error  equations 
presented  above.  The  agreement  obtained  is  excellent  except  in  the  single  case  of  the 
effect  of  camera  heading  angle  errors  on  store  attitude  (equation  (47)).  The  effects  of 
pitch  and  yaw  interact  markedly  here  but  even  so  the  simple  formula  presented  holds 
with  *50%  for  elevation,  heading  and  depression  angles  in  the  range  ±30°  and  this  is 
adequate  for  most  design  purposes. 

There  are  of  course  many  other  aspects  which  must  be  considered  in  arriving  at  an 
overall  specification  of  system  errors  such  as  the  problems  of  distortion  in  the  optics, 
lens  focus,  depth  of  field  and  resolution  and  vibration.  These  however,  are  beyond  the 
scope  of  this  Report. 
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5  COM  CL  US  IOWS  j 

i 

the  derivation  of  the  equations  relating  the  position  and  attitude  of  a  cylindrical 
object  to  the  sight-line  tangent  angles  Measured  in  two  vertical  planes  has  been 
presented  together  with  an  iterative  technique  for  their  solution.  The  convergence  of 
the  iterative  technique  has  been  illustrated  by  results  obtained  with  a  programme 
written  for  a  TI59  calculator. 

Simple  first  order  approximations  have  been  derived  for  the  dependence  of  attitude 
and  position  errors  on  both  measurement  errors  and  errors  in  setting  up  the  camera 
geometry.  It  is  shown  that  the  resultant  attitude  errors  can  be  very  many  times  greater 
than  the  measurement  errors  causing  them  and  that  the  technique  is  best  suited  to  long 
stores,  enabling  a  wide  separation  of  the  camera  planes,  and  for  relatively  short  range 
viewing  (R^/a  small) . 
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Appendix  A 

THE  EQUATION  OF  THE  STOBE  SURFACE 

Consider  the  point  F  on  the  surface  of  a  cylindrical  store  of  radius  a  with 
coordinates  Rad  in  the  coordinate  system  shown  in  Fig  2  and  having  its  origin  at  a 
point  0  on  the  axis  of  the  store.  Angle  a  represents  the  elevation  angle  of  point  P 
above  the  horizontal  XT  plane  while  &  is  the  azimuth  angle  measured  in  this  horizontal 
plane  relative  to  the  X  axis. 

The  elevation  and  azimuth  angles  defining  the  attitude  of  the  store  are  6  and  \ p 
respectively,  in  accordance  with  convention. 

Point  C  on  the  axis  of  the  store  is  the  centre  of  the  normal  circular  section 


whose  circumference 

contains  the  point 

P  . 

Thus 

PC 

«  a 

(A-l) 

OP 

-  R 

(A-2) 

and  by  Pythagoras: 

OC 

■  (*2  -  4  ■ 

(A-3) 

CD  and  PA 

are  perpendiculars  from  C  and  P  respectively,  to 

the  horizontal 

XY  plane,  CB  is 

the  perpendicular  from  C  to  PA  as  shown  in  Fig  2. 

Thus 

AB  -  CD  • 

(r2  -  a2)  sin  8 

(A-4) 

PA  - 

R  sin  a 

(A-5) 

and 

PB  -  PA  -  AB  - 

12  2\^ 

R  sin  a  -  l R  -  a  1  sin  8 

(A-6) 

In  triangle  DOA 

OD  - 

(r2  -  a2f  cos  8 

(A-7) 

OA  - 

R  cos  8 

(A-8) 

DOA  - 

B  -  ip  • 

(A-9) 

By  the  cosine  rule 

2  2 

AD  -  OD  + 

OA2  -  2  .OD  .OA  cos  (DOA)  . 

(A- 10) 

But 

AD  -  BC 

(A-l 1) 

and,  by  Pythagoras: 


2  2  2 
BC  +  PB4  -  a 


(A- 12) 
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hence 


|r  sin  a  -  (r2  -  a2f  sin  e|  +  (R2  -  a2)  cos20  +  R2  cos20 


A 


-  2r(r2  -  a2)  cos  a  cos  0  cos  (8  -  <l>) 


(R2  _  a2)  -  2R^R2  -  a2j  jsin  a  sin  0  ♦  cos  a  cos  0  cos(8  -  i|>)| 
and  for  R  *  a: 


2  2  21  . 

R  -  a  -  R  sm  a  sin  0  +  cos  a  cos  0  cos 


(8  -  *)| 


or,  rearranging,  we  have  for  Che  equation  of  the  surface  of  the  store: 

o  .2 


R 


I  -  {sin  a  sin  0  +  cos  a  cos  6  cos (8  -  <i)  { 


2  ‘ 


a 2  (A-J3) 

0  (A- 14) 


(A-15) 


(A- 16) 
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Appendix  B 
TOE  ELLIPSE 

The  familiar  equation  of  an  ellipse  in  Cartesian  coordinates  is: 

2  2 

+  *2  -  J  (B-I) 

P  <1 

where  p  and  q  are  the  magnitudes  of  the  semi-axes  in  the  x  and  y  directions 
respectively. 

Transforming  to  polar  coordinates  by  means  of  the  substitutions 

x  »  R  cos  4 
y  »  R  sin  4 


(B-2) 

(B-3) 


gives  the  polar  form  of  the  ellipse: 


2  2  2.2 

R  cos  4  ,  R  sin  4 

2  2 
P  q 


which  can  be  rearranged  tc*  give: 


(B-4) 


(B-5) 


If  we  now  rotate  this  coordinate  system  through  an  angle  ft  such  that  the  angle  a 
measured  in  the  new  system  is  given  by 


a  -  4  -  ft 


(B-6) 


we  have 


(B-7) 


This  is  seen  to  have  exactly  the  same  form  as  equation  (2)  of  the  main  text  which 
therefore  represents  an  ellipse. 

Comparing  coefficients  we  have 

a  •  p  (B-8) 

2 

K  -  1  -  Ey  (B-9) 

q 


ft  -  Y  . 


(B-IO) 
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Hence 


and  since 


1  *  - 2 

0  -  K) 

0  <  K  <  1 

>  a 


a 


0  -  K) 


(B-l 1) 


then  equation  (2)  therefore  represents  an  ellipse  with  minor  semi-axis  a  and  major 
semi-axis  a! (I  -  K)^  ,  rotated  through  an  angle  y  . 
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DEVELOPMENT  OF  THE  SUM  AND  DIFFERENCE  EQUATIONS 


C.  J  The  sum  equation 

From  the  sum  of  the  sight-line  tangent  angles  to  the  top  and  bottom  of  the  store 
given  in  equation  (9)  of  the  main  text  we  have: 


tan(a  j  +  a2)  »  tan^-|  -  cj  +  2y  ^ 


where  C  ■  tan 


(1  -  K)  cos  2(an  -  y)  +  K  -s- 


U  -  K)  sin  2  (aQ  -  y) 


expanding: 


tan(aj  +  a^) 


cot  C  *  tan  2y 
1  -  cot  C  tan  2y 


2 

(1  -  K) jsin  2(a Q  -  y)  cos  2y  +  cos  2(aQ  -  y)  sin  2y^  +  k2j  sin  2y 

_ R0 

l  )  a2 

(1  -  K)|cos  2(0^  -  y)  c°s  2y  -  sin  2(aQ  -  y)  sin  2y|  +  K— ^  cos  2y 


8in  2a0+T41[^sin2Y 


K  CL 

cos  2aQ  +  |  K  — j  sin  2y 
*0 


Hence  rearranging: 


siniOj  +  a2)  cos  2oq  -  cos(aj  +  a2>  sin  2aQ  R  fl2 


cos  (a  t  +  a2) 


-  j-  R  K  cos  2yjtan  2y  -  tan(aj  +  “2^} 
*0 

. (C— 5) 


2 

.in{(a,  ♦  «2>  -  2aQ}  -  %  ,in{2y  "  (“l  *  °2)} 

*0 


-  y  (a,  +a2)  -  ^  sin{2y  -  (a,  ♦  a2>}  . 
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C.2  The  difference  equation 

From  the  difference  between  the  sight-line  tangent  angles  given  in  equation  (10)  of 
the  main  text  we  have: 


2  S 

cos  (a,  -  a.)  “  -5 - t 

1  1  P  +  Q 


(C-8) 


,  ,  4  2 

(1  -  K)z  +  (2  -  Kr  %  -  2(2  -  K)  ( 1  -  K) 

_ ^0 _ % 

2  n2  5 

(1  -  K)  +  2K(1  -  K)  cos  2(aQ  -  o)  +  K =~r 


*0 


r: 


2  ,-2  _ 


which  yields  a  quadratic  in  a  /Rq : 


(C-9) 


4  2  (  > 

^^4(1  -  K)  +  K2  sin2(a2  -  a^J  -  2  (1  -  K)j2  +  K^cos2(a2  -  Oj)  cos  2(aQ  -  y)  -  iM 

«0  *0  '  ’ 

+  (1  -  K)2  sin2(a2  -  Uj)  =0  .  (C-10) 

Solving  the  quadratic,  choosing  the  correct  root  by  inspection,  and  after  some  tedious 
manipulation  we  obtain: 


2  +  K {cos  (a2  -  Oj)  cos  2(Oq  +  y)  -  1^  —  2  cos(a2  -  c^)  x 

x  |l  -  K  sin2(aQ  -  y)^2  -  K  +  K  cos2(oq  -  y)  cos2(aQ  -  “j)^' 


4 


4  +  j~-j^  sin2(a2  -  Oj) 


(C-ll) 
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Appendix  D 
CAMERA  (EOMETRY 

Fig  4  shows  the  basic  geometry  and  relative  positions  of  the  store  centres  in  the 
two  camera  planes.  The  camera  separation  is  0  and  the  separation  of  the  store  centres 
measured  along  the  axis  of  the  store  is  L  . 

Resolving  into  the  horizontal  and  vertical  planes  as  shown  in  Fig  4  we  have,  from 
the  horizontal  plane: 


and 


L  cos  8  cos  ip  =  D  +  Rq  cos  <Jq  cos  ipQ  -  Rq  cos  a q  cos 

L  cos  8  sin  ip  =  Rq  cos  sin  ip^  -  Rq  cos  a q  sin  i|Jq 


(D-l) 

(D-2) 


and  from  the  vertical  plane: 


L  sin  0  = 


Ro  sin  a0  -  Ro  sin  “o 


Eliminating  L  from  (D-l)  and  (D-2)  gives: 


Rq  cos  aQ  sin  ip0  -  Rq  cos  ct’Q  sin  ^ 
tan  *  D  +  Rq  cos  aQ  cos  (Pq  -  Rq  cos  a'Q  cos 


Similarly  from  (D-2)  and  (D-3) 


tan  6  =  sin  p 


Ro  sin  °0  '  Ro  sin  “o 


|Rq  cos  aQ  sin  ip 0  “  Rq  cos  “q  sin 


(D-3) 


(D-4) 


(D-5) 


where  the  sin  iJj  term  on  the  right  hand  side  may  be  eliminated  by  substitution  using 
equation  (D-4)  if  required. 

Equations  (D-4)  and  (D-5)  give  the  store  attitude  in  terms  of  the  camera  geometry 
and  positions  of  the  store  centres  and  they  may  be  rearranged  to  give  the  relationship 
between  the  positions  of  the  store  centres  for  a  known  attitude. 

Rearranging  (D-5)  we  have: 

(  sin  8  cos  cos  o  +  cos  8  sin  1 1>  sin  a 
%  "  *0)  sin  8  cos  i|»q  cos  c*q  +  cos  8  sin  ip  sin  clq 

and  from  (D-4) 

D  sin  ip  +  R  cos  Oq(cos  «Pq  sin  ifr  -  sin  ^  cos  (p> 

-  R'  cos  <*q(cos  i|/q  sin  p  -  sin  cos  p)  .  (D-7) 


28 


Appendix  D 


This  simplifies  to  give: 


cos  Oq  s  in  (if'.,,  -  ifr)  -  D  sin 
cos  Qq  sin(ifiQ  -  ijii  ~ 


Equating  the  right  hand  sides  of  equations  (D-6)  and  (D-8)  and  collecting 
in  cXq  yields: 


tan  c*q 


Rq  sin  6  cos  aQ  sin(ifiQ  -  \j/g)  +  cos  0  sin  Oq  sin(#g  -  t)|  •  D 


3 


cos  8<R_.  cos  a  sin^.. 


*) 


D  sin 


♦} 


(D-8) 

together  terms 

sin  6  sin 

.  (D-9) 
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Appendix  E 

DETERMINATION  OP  BANK  ANGLE 

Consider  a  line  painted  along  the  port  side  of  the  cylindrical  store  (the  side 
assumed  to  be  nearest  to  the  cameras),  so  that  in  straight  and  level  flight  the  line  lies 
in  the  horizontal  plane  containing  the  centre  line  of  the  store.  The  problem  is  to 
determine  the  sight-line  angle  of  the  painted  mark,  in  either  camera  plane,  when  the 
store  banks  through  an  angle  $  defined  in  the  conventional  sense  (port  wing  high), 
and  is  otherwise  of  known  attitude  and  position. 

A  procedure  may  be  followed  which  is  similar  to  that  adopted  in  determining  the 
sight-line  tangent  angles  to  the  top  and  bottom  edges  of  the  store.  The  general  equation 
of  the  painted  line  in  three-dimensional  space  is  found  first  and  the  intersection  of 
this  line  with  the  camera  viewing  plane  then  gives  the  coordinates  of  the  observed 
point  where  sight-line  angle  is  required.  The  foregoing  is  carried  out  in  a  polar 
coordinate  system  centred  on  the  store  centre  line  and  a  simple  transformation  of  axes 
to  a  new  origin  at  the  camera  position  enables  the  sight-line  angle  to  be  found  in  the 
coordinate  system  used  for  experimental  measurements . 

The  geometry  of  the  situation  is  shown  in  Fig  5  where  P  is  a  point  on  the  line  at 
coordinates  Ra8  ,  a  and  3  being  the  elevation  and  azimuth  of  the  point  P 
respectively.  C  is  the  centre  of  the  normal  section  having  P  on  its  circumference  and 
CH  is  the  horizontal  from  C  in  the  same  normal  section  of  the  store. 

PA  and  CD  are  the  perpendiculous  from  P  and  C  respectively  on  to  the 
horizontal  XY  plane  containing  the  origin  0  ,  which,  as  indicated  earlier  lies  on  the 
centre  line  of  the  store  (OC) .  Also,  PQ  is  the  perpendicular  from  P  to  CH  , 
giving  APQ  »  0  ,  AE  is  the  perpendicular  from  A  to  OD  ,  QJ  is  mutually 
perpendicular  to  PQ  and  CH  with  J  lying  on  PA  ,  and  JK  is  the  perpendicular  from 
J  to  OC  . 

Now,  since  QJ  and  CK  are  both  mutually  perpendicular  to  the  normal  section 
PHC  ,  QJKC  is  a  rectangle  and  therefore  JK  is  equal  to  QC  and  also  horizontal.  It 
follows  that  JKAE  is  also  a  rectangle.  We  have: 


PA  -  R  sin  a  -  PJ  +  AJ  (E-l) 

-  PQ  sec  0  +  KE  (E-2) 

■  a  sin  $  sec  0  +  OE  tan  0  (E-3) 

■  sec  0{a  sin  0  +  OA  cos(<>  -  8)  sin  of  (E-4) 

therefore  R  sin  a  -  sec  ©{a  sin  9  ♦  R  cos  a  cos  (ifr  -  0)  sin  e|  .  (E-5) 


Also: 

hence 

> 

) 


AE  -  OA  sin(4>  -  0)  -  JK  -  QC 
R  cos  a  «in(t|>  -  0)  ■  a  cos  9  . 


(E-6) 

(E-7) 


L 
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Substituting  from  (E-7)  for  R  cos  a  into  (E-5)  gives  two  independent  equations 
in  R  and  a  which  nay  be  solved  for  these  parameters  if  required,  however  this  is  not 
necessary  for  present  purposes 


R  sin 


afsin  $  sin(^  -  0)  +  sin  0  cos  4  cos(iP  -  B)> 

in  a  -  ' -  ■  —  ■  ■ — r—n - ct- - J 

cos  0  sin(ip  -  B ) 


(E-8) 


R  cos  a 


a  cos  A 
sin<4>  -  6)  * 


(E-9) 


As  in  determining  the  equation  of  the  store  section  in  the  camera  plane 
(section  2.2)  we  now  assume  that  the  camera  plane  passes  through  the  origin  with  azimuth 
,  enabling  the  coordinates  of  the  mark  observed  in  the  camera  plane  to  be  determined 
by  substituting  B  m  >l>g  +  n  in  equations  (E-8)  and  (E-9)  above  for  the  port  side, 

6  *  ♦  o  giving  the  starboard  side.  If  we  also  change  the  coordinate  system  as  before, 
giving  a  new  origin  at  the  camera  position  and  specifying  the  new  coordinates  of  point  P 
as  range  and  depression  angle  ,  we  have  for  the  port  side 


Rj  cos  ■  Rq  cos  Oq  -  R  cos  a  (E-10) 

R2  sin  a3  -  Rq  sin  oQ  -  R  sin  a  (E— 11) 

where  Rq  and  aQ  are  the  range  and  depression  angle  of  the  store  centre  in  the  new 
coordinate  system. 

Hence: 

R0  sin  aQ  -  R  sin  a 

tan  a  -  - - - - - - 

3  Rq  cos  aQ  -  R  cos  a 

and  substituting  from  equations  (E-8)  and  (E-9) 


tan  a. 


U  . 

—  sin  QqC08  0  sin(i|iQ  -  ^)  -  sin  $  sin(ip0  -<(»)+  sin  0  cos  $  cos(tt>0  ♦  ip) 
cos  cos  Uq  sin(i|>Q  -  <p)  -  cos 


(E— 12) 


Note  that  Rq  has  been  normalized  by  the  store  radius  a  and  the  angle  (<p  -  iPq) 
changed  to  (tpQ  -  <p)  to  agree  with  previous  form. 

Equation  (E-12)  may  be  rearranged  to  give  the  bank  angle  $  as  a  function  of  the 
sight-line  angle  a^  by  cross  multiplying  and  then  collecting  together  terms  in  sin  ♦ 
and  cos  4  ,  thus: 

^cos  0  tan  a3  +  sin  0  cos(ip0  -  <p)  J  cos  ♦  -  sin(<pQ  -  ip)  sin  $ 


•  —  cos  0  sin  ( ip. 


-  ip){ 


tan  a,  cos  a.  -  sin  a. 


(1-13) 
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Table  2 

CONVERGENCE  OF  ITERATIVE  PROGRAMS  FOR  30°  PITCH  AND  YAM 


Input  data 

D 

a 

1.25 

°1 

•  9.814760222 

*0  ■ 

o 

o 

o\ 

“2 

-  38.9993655 

*Q  ’ 

vo 

o 

0 

-  21.01570820 

°2 

-  51.88559048 

Results 

Parameter 

number 

0 

* 

“0 

Jo 

a 

“0 

Jo 

a 

0 

23.178612 

25.349156 

24.407063 

3.969209 

36.450649 

3.757364 

1 

28.105158 

27.478669 

23.857871 

4.646414 

36. 174493 

4.458827 

2 

29.074128 

29.269900 

23.953318 

4.937227 

36.314031 

4.730071 

3 

29.704085 

29.595812 

23.935244 

5.051800 

36.294045 

4.847901 

A 

29.847944 

29.877479 

23.94  51 

5.099939 

36.309429 

4.892836 

5 

29.950265 

29.932039 

23.944922 

5. 119082 

36.306468 

4.912509 

6 

29.974322 

29.979157 

23.946859 

5.127210 

36.308892 

4.920093 

7 

29.991542 

29.988483 

23.946448 

5.130457 

36.308409 

4.923427 

8 

29.995643 

29.996446 

23.946770 

5.131837 

36.308811 

4.924716 

9 

29.998559 

29.998045 

23.946702 

5. 132389 

36.308731 

4.925282 

10 

29.999260 

29.999394 

23.946756 

5. 132624 

36.308799 

4.925502 

00 

30.0 

30.0 

23.946753 

5.132785 

36.308797 

4.925627 

JL  cT 
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LIST  OF  SYMBOLS 

store  radius 
camera  separation 
store  vertical  drop  height 
parameter  defined  by  equation  (2) 

as  K  but  referred  to  aft  camera  plane  (see  equation  (28)) 
semi-axis  of  ellipse 

variable  expression  defined  by  equation  (6) 
semi-axis  of  ellipse 

variable  expression  defined  by  equation  (6) 
radius  from  origin  on  store  centre  line 
radius  from  origin  at  camera  position 
range  of  store  centre  from  camera  position 
as  Rq  but  referred  to  aft  camera  plane 
range  of  mark  on  store  side  from  camera  position 
variable  expression  defined  by  equation  (6) 
general  coefficient  of  cos  4  in  equation  (E-14) 
general  coefficient  of  sin  $  in  equation  (E-14) 
general  elevation  angle  from  origin  at  store  centre 
general  depression  angle  from  origin  at  camera 
depression  angle  of  store  centre  from  camera  position 
as  aQ  but  in  aft  camera  plane 

depression  angle  of  sight-line  tangent  to  top  of  store  from  camera  position 
as  a(  but  in  aft  camera  plane 

depression  angle  of  sight-line  tangent  to  bottom  of  store  from  camera 
position 

as  0^  but  in  aft  camera  plane 

depression  angle  of  mark  on  camera  side  from  camera  position 

general  azimuth  angle  from  origin  at  store  centre  measured  relative  to 
direction  of  line  joining  the  two  cameras 

parameter  defined  by  equation  (2) 

as  y  but  referred  to  aft  camera  plane  (see  equation  (28)) 
increment  in  range  RQ 

increment  in  range 


LIST  OF  SYMBOLS  (concluded) 


increment  in  depression  angle 
increment  in  depression  angle  a ^ 
increment  in  depression  angle  a( 
increment  in  depression  angle  aj 
increment  in  depression  angle 
increment  in  depression  angle  n^ 
increment  in  store  elevation  angle  6 
increment  in  camera  elevation  angle 
increment  in  store  heading  angle  \p 
increment  in  camera  heading  angle  i|» 
increment  in  aft  camera  heading  angle  ^ 
mean  increment  in  store  range 
mean  increment  in  store  depression  angle 
angle  defined  by  equation  (E— 14) 
store  elevation  angle 
angle  defined  by  equation  (7) 
store  bank  angle 

polar  coordinate  defined  by  equations  (B-2)  and  (B-3) 

store  heading  angle 

forward  camera  heading  angle 

aft  camera  heading  angle 

angle  used  in  equation  (B-6) 
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